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Abstract 

We present microscopic calculations of neutrino propagation in hot neutron 
matter above nuclear density within the framework of the Random Phase Ap- 
proximation. Calculations are performed for non-degenerate neutrinos using 
various Skyrme effective interactions. We find that for densities just above 
nuclear density, spin zero sound is present at zero temperature for all Skyrme 
forces considered. However it disappears rapidly with increasing temperature 
due to a strong Landau damping. As a result the mean-free path is given, 
to a good approximation, by the mean field value. Because of the renormal- 
ization of the bare mass in the mean field, the medium is more transparent 
as compared to the free case. We find, in contrast, that at several times nu- 
clear density, a new type of behavior sets in due to the vicinity of a magnetic 
instability. It produces a strong reduction of the mean free path. The corre- 
sponding transition density however occurs in a region where inputs from more 
realistic calculations are necessary for the construction of a reliable Skyrme 
type parametrization. 

PACS numbers: 11.30. Qc, 05.70.Fh, 12.38.Mh, 25.75.+r 
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I. INTRODUCTION 



The gravitational collapse of massive stars at the end of their thermonuclear burning is 
believed to produce a core of hot dense matter in which most of the initial gravitational 
binding energy is stored into neutrinos. Numerical simulations of the subsequent evolution, 
leading to the formation of a neutron star, require the knowledge of the equation of state of 
hot dense matter as well as a reliable description of neutrino transport phenomena. Relevant 
quantities are the specific heat of the medium and the neutrino mean free path. 

The mean free path of a neutrino due to scattering inside neutron matter at temperature 
T is proportional to the optical potential. It can be expressed in the case of non-degenerate 
neutrinos as ill 
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where Gf is the Fermi constant, qa the axial coupling constant, kj and kj are the initial 
and final neutrino momenta, q is the transferred momentum kj — k/, u is the transferred 
energy |kj| — |k/|, and cos 9 = kj • k/. The functions 5^(a;,q, T) represent the dynamical 
structure factors in the spin symmetric (5=0) or spin antisymmetric (5=1) channels. They 
are given by the imaginary part of the respective response function x (^j Q) an d contain 
the relevant information on the medium. The vector (axial) part of the neutral current 
gives rise to density (spin-density) fluctuations, corresponding to the 5=0 (5=1) channel. 
Although we shall consider only nondegenerate neutrinos, we note that the previous formula 
should be corrected for additional Pauli blocking factors in the degenerate case. 

Several authors have evaluated the neutrino mean free path using various approximation 
schemes and various models of the trapping environment. In their early calculations Tubbs 
and Schramm omitted the effect of nucleon-nucleon interactions, whose importance was 
later on emphasized by Sawyer || showing that for standard equations of state, the neutrino 
opacity due to coherent scattering off density fluctuations can be significantly larger than 
in a neutron Fermi gas. Iwamoto and Pethick investigated the effect of nucleon-nucleon 
interaction within the framework of Landau theory of Fermi liquids, considering both density 
and spin-density fluctuations according to Eq. ([]]). The effective interaction was taken into 
account by the monopolar Landau parameters F and Go derived earlier by Backman et al. 
||]. The main conclusion of their work is that for degenerate neutrinos, there is a reduction of 
in-medium scattering cross sections by a factor of 2-3 due to nucleon-nucleon interactions, 
and a subsequent increase of the neutrino mean free path. Calculations have also been 
performed using the ring approximation in which the response function X^O^o) °f 

the interacting system is expressed in terms of the bare response function Xa{u, q) as 

X™M= , v lt\ q \ y (2) 

l - y( 6 )(q)xo(^,q) 

where V^ 5 ^(q) is the interaction in momentum space in spin channel 5. 

The case of non-degenerate neutrinos was examined by Haensel and Jerzak ||, who 
computed both the scattering and the absorption rates within the framework of Landau 
theory. They used a quasiparticle interaction derived from the Reid soft core nucleon-nucleon 
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potential. Their results indicate that the higher the density of the medium, the larger the 
effect of the interaction. More recently, Reddy et al || considered degenerate neutrinos 
immersed in an environment of nucleons interacting via a Skyrme two-body interaction out of 
which the monopole Landau parameters F and Go are constructed. These authors consider 
the Skyrme parametrizations SGII and SkM* and simulate the momentum dependence of 
the force through the direct matrix element of the meson exchange potential. 

In this paper we extend preliminary results presented in Ref. and report systematic, 
fully self-consistent calculations based on effective interactions of the Skyrme type, including 
some of the most recent ones. Such interactions have been very successful at describing 
properties of finite nuclei and nuclear as well as neutron matter. An attractive feature of 
Skyrme interactions is that response functions can be calculated in closed form IPHTT 



in 



the Random Phase Approximation (RPA), which generalizes the simple ring approximation 
(@) based on Hartree-like dynamics, to the more realistic case of Hartree-Fock mean field 
dynamics. An advantage of such a framework is that it also avoids splitting structure factors 
into single particle and collective contributions based on different approximation schemes. 
It further avoids postulating a factorization of the temperature dependence of the structure 
factors, which are in conflict with the fact that zero sound disappears at high temperature 
I"I2| ] . Moreover, sum rules hold in the case of fully self- consistent calculations, thus allowing 



the possibility of consistency checks. 

In the following we focus on the evaluation of the mean free path of neutrinos interact- 
ing with neutrons through neutral currents only. This coupling is important mainly after 
deleptonization has taken place and it affects the late behaviour of the neutrino burst, where 
neutrino diffusion is driven by the temperature gradient. Indeed, in the earlier stages, just 
after the formation of the protoneutron star, charged current reactions leading to neutrino 
absorption are dominant, due to the fact that the high temperatures involved prevent ef- 
ficient inhibition of charged reactions by the energy-momentum conservation constraints. 
We consider neutron matter above nuclear density po=0.17 nucleons/fm 3 . The case of high 
densities is of special interest because medium effects are expected to play an important 
role. In contrast, below nuclear density propagation is dominated by coherent scattering off 
nuclei, for which reliable estimates of the opacities are available. 

This paper is organized as follows. In Sec. II we briefly summarize the formalism for the 
computation of the response of the neutron liquid and analyze the main results in Sec. III. 
Our conclusions and perspectives are presented in Sec. IV. 

II. RESPONSE FUNCTION OF HOT NEUTRON MATTER 



In the present work we use the techniques developed in Refs. [|10|,|TT1 to obtain RPA 



susceptibilities at finite temperature, in the particular case of Skyrme-type effective interac- 
tions. The RPA dynamical susceptibility x^Kq^^T) in spin channel S of a Fermi system 
in thermal equilibrium at temperature T is given by the the polarizability to an infinitesimal 
external field 

Vext = eO exp(iut — zq.x). (3) 

In this equation uj contains a small imaginary part to ensure an adiabatic switching of the 
external field starting at time t = — oo when the medium is in its equilibrium state at 
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temperature T. The operator O is equal to the unit matrix in spin space for spin zero (S=0) 
and to the third Pauli matrix a z for S=l. At late enough times (and small enough e) the 
expectation value of the operator O has the same space and time dependence as the external 
field. The dynamical susceptibility is defined by 
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lim 



(O) 



e^o e exp(iut — zq.x) ' 
and it can be shown to be given by the algebraic expression 

AS) 



r 



Xo + G V ph G , 



(4) 



(5) 



where Xo is the response of the free system, G is the free particle-hole (ph) propagator at 
temperature T and V p ^ is the RPA induced interaction, related to the effective ph one V^f 
by the integral equation 



V p f + V£> G V ph 



r(S) 



(6) 



In Ref. ||, it has been shown that for the general class of ph interactions of the following 
form containing a monopole and a dipole term 



Vl%i, g 2 ; q) = Wl S \q) + Wi S ) (g)(q 1 - q 2 )^ 



(7) 



to which the Skyrme ones belong, the solution of the integral equation can be algebraically 
obtained in a closed form, in terms of generalized Lindhard functions. This result was derived 
for the zero temperature case. However the generalization to a thermally excited liquid is 
straightforward. The effect of temperature appears implicitly in the RPA equations through 
the Fermi-Dirac occupation numbers contained in the free ph propagator. The algebraic 
structure of both the RPA equation and its solution is preserved [10,11] and the dynamical 
susceptibility in channel (S) has the form 



X {S \q^,T) 



Xo(q,uj,T) 



DW{q,u,,T)' 

where the factor 2 takes into account the spin degeneracy and 
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(9) 



In this expression kp is the neutron Fermi momentum while m* denotes the nucleon effective 
mass in the mean field. 
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Since in pure neutron matter the isospin exchange operator P r reduces identically to the 
unit matrix, only particular combinations of the Skyrme interaction parameters ij and Xj, 
namely 



s 2 = t 2 (l + x 2 ), 



Si = *i(l -Xi) 

S3 = £ 3 (1 - x 3 ) 



(10) 



are required. Thus, one obtains the following expressions for the neutron effective mass 



h 2 



h 



2m* 2m 



o 1" o( Sl + 3s s)p, 



(11) 



(S) 

where p is the density of neutron matter, and for the coefficients W[ 2 of the effective ph 
interaction 
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with 7 being the exponent of the density dependent term in the Skyrme force 
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(12) 
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In the special case of the Skyrme interaction SIII we have considered it as a three-body 
force rather than a density dependent one, i.e. we have used a value S3 = in the previous 
equations. 

The functions X2i(u,q) (i—0, 1, 2) are generalized susceptibilities per unit volume Q, 
defined as [01 



X2i 



n * 2 



n(k) 
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(k + q) 



(1 -n(k + q)) 
1 \ 



with 



u; — o;(k, q) + ir/ u; + u;(k, q) + ir) 



n(k) 



l + exp(e(k) -n)/T 



(14) 



(15) 



the Fermi-Dirac occupation probability, p being the chemical potential at the given temper- 
ature. Here 



w(k,q) = e(k + q) -e(k), 



(16) 
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with 

e(k)= h 2 k 2 /2m\ (17) 

the single-particle energy. The dynamical structure factor is obtained, at both positive and 
negative energies u, from the detailed balance relationship 

S( ] ^u,T) = _ e _^ /T . (18) 

It is worthwhile stressing that opposite to the ring approximation defined by equation (|J), 
which involves only the Linhard function xoi new terms involving the generalized suscepti- 
bility X2 an d Xi arise in the full RPA approximation. 

Another important point about our self- consistent approach is that it preserves sum 
rules. In particular one has the following energy weighted sum rule in the spin one case 

r+°° T) 2 1 

j ^ S^(q,u 7 T)udw = {^ + g(*i " s 2 )p}pq 2 . (19) 

In the spin zero case the right hand side is simply h 2 q 2 p/2m. 

In the following discussions regarding the appearance of zero sound modes it will turn 
out to be convenient to have the explicit expressions of the real and imaginary parts of the 
bare response function. Its imaginary part is given by 



m 



2 



kT , 1 + e ^ A +^ 



Zm X o(u>, q) = -^Y^ZJ log i - e , (A _, /2) . (20) 



In this equation (3= 1/kT and the quantity A is given by 



2 2 

muj q A . . 

lq l 8m 2 

The expression of the real part is a little more difficult to construct. In the zero temperature 
case it reads Jl7j 

$te X o(u, q, k F , T = 0) = + v {x + ) - y{xJ)}. (22) 

In this equation the quantities x± are defined by 

mu q 

X± = nk^q ± 2k- F ^ (23) 

and the function ip reads 

<p(x) = (l-x 2 )log|^|. (24) 
x + 1 

For a non zero temperature the real part of the response function is obtained by perform- 
ing an average of the zero temperature results obtained for various values k of the Fermi 
momentum with a weight factor equal to the occupation number n(k): 

3?e%o(^; Q, kp, T) — — / dtexo(^,q,kp = k,T = 0)dn(k). (25) 

Jo 
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III. RESULTS 



We have considered two standard Skyrme parametrizations, namely SIII (T3| and SkM* 
T4| , as well as parametrizations SLy230a and SLy230b recently adjusted |15] to reproduce 



the latest equation of state for pure neutron matter calculated by Wiringa et al [IB] within 
the Fermi Hypernetted Chain scheme, using the Urbana V14 two-body interaction plus 
the Urbana three-body interaction. The energies per particle given by parametrizations 
SLy230a,b and by the realistic force are very close to each other up to densities of about 1 
nucleon/fm 3 . In this section we analyze in detail the evolution of the mean free path with 
neutron density, its dependence upon temperature and neutrino momentum and the effect 
of residual interactions, as follows. 



A. Evolution of the mean-free path with density 

In Fig. 1 we show the evolution of the mean free path with density for an incoming neu- 
trino momentum hi— 5 MeV/c and a temperature T=5 MeV. Results are displayed for the 
four effective interactions mentioned above. It can be noticed that for interactions SLy230a 
and SLy230b the mean free path exhibits a slow variation with density. This behavior is 
similar to the one obtained by Haensel and Jerzak 0] who find mean free paths A= 2, 1.5 
and 1km for densities 0.2, 0.36 and 0.6 fm -3 respectively, for fcjC = T = 5MeV. In contrast 
the results obtained with interactions SIII and SkM* show a substantially different behavior, 
both exhibiting a rapid monotonic decrease with density. This decay is the consequence of 
a pronounced increase in the magnetic strength, that reflects in the vicinity of an instability 
in the channel S=l. The magnetic instability can be located by examining the Landau 
parameters of neutron matter, in particular, the Landau coefficient Go related to the spin 
asymmetry coefficient a a in neutron matter through the relation 



h 2 k 2 F [ 

6m* 



l + G ). (26) 



In terms of the force parameters, the asymmetry coefficient is given by the following 
expression 

a - = ~ s o + — P 7 h + — [S2 - si r. 27 
6m* 6 4 24 

where r = 3/5 p k F is the kinetic energy density. For interaction SIII the last term dominates 
at high values of the neutron density and as a result the Landau parameter takes the value -1 
at a critical density around 0.317 fm~ 3 . This critical density p c is (in fm -3 ) 0.1967 for SkM*, 
0.5411 for SLy230a and 0.5928 for SLy230b. These are precisely the abscissae at which the 
neutrino mean free path drops to a very small value in Fig. 1. Actually, the mean- free path 
at the minimum is not exactly zero, since even though the monopole contribution vanishes 
there remains a small contribution from the dipole term. 
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B. Spin zero-sound with Skyrme forces 



To illustrate this point, in Fig. 2 we display the dynamical structure factors computed 
with interaction SIII, for spin channels S=0 and 1 (left and right parts of the figure, re- 
spectively), as functions of energy lu, for transferred momenta q equal to 5 and 100 MeV/c 
(dashed and solid lines, respectively) and a temperature of 5 MeV. The upper and lower 
parts in this figure correspond to densities 0.1 fm~ 3 and p c . The zero sound mode is visible 
in the S=l, p= 0.1 fm~ 3 case although it is smeared by Landau damping at a temperature 
T=5 MeV. We also clearly appreciate the significant increase and change of scale in S^^u;, q) 
as the density reaches the transition point, which contrasts with the rather moderate vari- 
ation undergone by S 1 - ^ with density. The numerical data indicate a scale factor of three 
orders of magnitude in as the density reaches p c . For S=0 one also observes a peak 
which corresponds to the maximum in the imaginary part of the bare response. This peak is 
reduced and shifted towards lower energies in the RPA response because of the contribution 
of the real part. 

In order to have a qualitative criterion for the existence of a zero sound mode it is 
convenient to calculate the response in the monopole approximation. In this approximation 
one evaluates the particle -hole interaction strength at the Fermi surface and one ignores its 
angular dependence i.e. one sets q\= q\= k 2 F and qi.q2=0 in Eq. (|?p. One thus uses Eq. 
(|) replacing V^{q) by W 1 + 2W 2 k 2 F i.e. 

Vi S=1) = so - \s,p~< + \{s 2 - s x )k 2 F - hsi + s 2 )q 2 , (28) 
2 4 

in the spin one case, and by 

Vf =0) = s + (7 + l)(7 + 2) g3p7 + 1 + 3sa)4 + _ 3s2)g2) (2Q) 

in the spin zero case. From the expression for the structure function 

q , rp^_ _! 25m Xo(^,q) 1 /o n x 

tt(1-V^ >$texo) 2 + {Vo '^mxo) 2i - e 

we see that, for a given value of the momentum q, zero sound will be present whenever there 
exists a solution ojr of the equation 

l-V (S) 3?e X oK,q) = 0, (31) 

in the energy domain where the imaginary part of the bare response function is small (or 
zero). At zero temperature this domain corresponds to (see equation(pO|)) 



h 2 a 

uJ R >u + = —±(2k F + q). (32) 
2m 

At this point the real part of the response function is given by equation(|22"D i.e. 

3?e XoM = -^L{(2k F + q) log(^±^) - 2k F }. (33) 
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Since the real part of the response function decreases with energy in the domain of interest 
we find that the condition for zero sound at zero temperature is 

V {S) ^e Xo ^ + ,q)>l. (34) 

For interaction SIII one finds that for a momentum q of the order of 5 MeV/c and a density 
of 0.1 neutron per fm 3 one has 

V Q {S=1) = 112.9MeV x fm 3 and V Q {S=0) = -560.8MeV x fm 3 (35) 

and 

Xo(^+, q) = 6.2 x lO^MeV" 1 x fm" 3 . (36) 

These values show that the residual particle-hole interaction in the S=0 channel has the 
wrong sign to produce a zero sound mode. In contrast in the S=l case it is repulsive and 
strong enough to give a zero sound mode. It is however not as strong as in the case of liquid 
helium-3 [|T7]] where zero sound occurs far in the quadratically decreasing tail of the real part 
of the response function. Here in contrast, the energy of the zero sound mode is close to 
the maximum of the real part of the response function, uj max = qkpjm. Note that this value 
belongs to the integration domain of equation (|1|), namely 

\u\<q<2E v -u u<E u , (37) 

thus allowing, in principle, the zero-sound mode to contribute to the neutrino mean-free path. 
It is worthwhile noting that for a strong particle-hole interaction the magnetic strength would 
be driven outside the integration domain which would significantly increase the mean-free 
path. 



C. Energy density near the magnetic transition 

At the critical density the symmetric Hartree-Fock ground state becomes unstable. This 
can be visualized looking at the expression of the energy density of spin polarized neutron 
matter: 

h 2 1 1 

U = 2^* ^ + T ^ + <yS ° + 6 S3P ^ P]Pi + 8^~ Sl + S2 ^ PT ~ Pi ^ T] ~ n " > ' ^ 

from which one can construct the energy per neutron as a function of the spin polarization 
parameter 

x = P ^- P ±. (39) 

The minimum energy per particle occurs at x=0 for densities smaller than the critical one 
and at nonvanishing asymmetries for larger densities. It actually reaches rapidly the value 
x=l which corresponds to a fully polarized medium. 

In the vicinity of the magnetic transition, the response function in the S=l channel 
behaves, according to equation (01), as the inverse of the spin asymmetry coefficient a a in 
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the limit oo=0 and q=0. This can be checked by noting that in the limit of a small momentum 
transfer one has 

Ke Xo (uJ, q = 0, k F , T = 0) = -^f , (40) 

and a negligible value of the imaginary part of the bare response function. The RPA response 
thus becomes [0 



X (S=1 V = 0,q=0) = --^. (41) 

Z,(X a 

This formula accounts for the presence of the peak which develops near zero energy in figure 
2 when the density becomes close to its critical value. It also shows that in the limit of 
a vanishing momentum transfer the mean-free path would exactly vanish at the transition 
point. This result is not in conflict with the sum rule mentionned previously (see eq. (|T9|)) 
because a peak at zero energy gives no contribution to the energy weighted sum. We have 
not attempted to push our RPA calculations beyond the critical point, since more work is 
needed to produce a reliable Skyrme type parametrization in this region. 



D. Momentum and temperature dependence 

In Fig. 3, the mean free paths obtained with the interaction SLy230b are displayed as 
functions of the density for six combinations of hi and T. It is worthwhile noticing that 
the effect of increasing momentum and/or temperature is similar and mostly consists of an 
overall reduction of the scale. Scaling properties were already discussed by Haensel and 
Jerzak |J. As a result of their approximation scheme these authors found that scaling the 
neutrino energy and the temperature by a factor a produces a scaling of the mean free path 
by a factor 1/a 3 . Comparing the curves fcjC = T— 5 MeV and fc^c = T— 10 MeV we see 
that this scaling law holds to a reasonable accuracy in our RPA calculations, although the 
curvature of the uppermost line is somewhat too high. 



E. Effect of residual interactions 

An important issue regarding neutrino propagation is the effect of the residual neutron- 
neutron interaction. Various answers to this question are found in the literature. Indeed 
while Iwamoto and Pethick conclude that Fermi liquid effects increase the mean free path 
to some extent, Haensel and Jerzak |5| find that near normal density, the overall effect of the 
nucleon -nucleon interaction on the mean free path is rather weak. In contrast, the recent 
paper by Reddy et. al. || reports that many-body effects suggest considerable reductions 
in the opacities compared to the free gas estimates. 

Our results are displayed for interactions SIII and SLy230b in the cases T = fcjC = 5 MeV 
and T = kiC = 10 MeV in Fig. 4. This figure shows the relative opacities X~ x / Xj ree (m*) 
where Aj r ! ee (m*) corresponds to a gas of free neutrons with effective mass m*, rather than m, 
i.e. the scale is furnished by the opacity in the mean field. In the vicinity of the saturation 
density of nuclear matter it can be seen that this relative opacity is close to unity. This means 
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that for such densities, the major role of the interaction is to dress the bare mass giving 
rise to the mean field density dependent effective one m*. If the reference for the relative 
opacity is however chosen as A/ ree (m), then one finds that interactions do make the medium 
more transparent. Indeed at normal density one has m*/m= 0.695 for interaction SLy230b 
and the smaller the mass, the smaller the opacity. In particular, at zero temperature 



Xf ree 2AtC 3 \m* + kp / 



(42) 



Since the effective mass becomes small at high density, large reductions can occur. Never- 
theless, at large enough densities the dominant effect in Fig. 4 is the reduction of the mean 
free due to the vicinity of the magnetic instability. 

The reason for having RPA values of the mean-free path very close to the mean field 
ones (at least in the region just above nuclear density), in spite of the presence of a zero 
sound mode at zero temperature, is that this mode is rapidly smeared as temperature 
increases because of the strong Landau damping, clearly visible on Fig. 2. Moreover when 
temperature rises there is a strong reduction of the real part of the response function, which 
may correspond to a factor 2-3 for a temperature of the order of 5 MeV (see e.g. figure 1 
of reference |18|]). As a result one finds that the term Vq S ~ 1 ^x m the expression of the RPA 
response function is substantially smaller than unity, which implies that RPA and mean field 
response functions are nearly identical in the density temperature range considered here. 



IV. DISCUSSION 

In the present work we have calculated neutrino propagation in neutron matter within a 
microscopic and consistent framework able to describe both the nuclear equation of state and 
dynamical structure factors. The approach adopted here is based on Skyrme type effective 
interactions and linear response (RPA) theory. As already mentioned, an advantage of 
the present procedure is the treatment of single particle and collective contributions to the 
structure factors on an equal footing at any temperature, making unnecessary to postulate 
a factorization of the temperature dependence of the structure factors. Such factorizations, 
which are in conflict with the fact that spin zero sound disappears at temperatures of a few 



MeV [|T2| , overestimate the contribution of collective modes. Our calculation indicates that 
magnetic zero sound modes appear at zero temperature for all Skyrme parametrizations here 
considered, even at the largest momenta (i.e., q = 100 MeV/c) up to density p c ; however, 
at finite temperatures they become subject to Landau damping on the ph continuum, as 
shown in Fig. 2 and discussed in section III. 

In the domain where the RPA is reliable (stable symmetric ground state in the mean 
field picture) we have found that the most important effects of the interaction are the renor- 
malization of the mass of the neutrons into their effective mass m* in the response function. 
This makes the medium more transparent, thus giving rise to shorter deleptonization time 
scales 0. 

We have found that at high density instabilities occur, which signal a breakdown of the 
Skyrme like parametrization. The possible existence of a magnetic transition in neutron 



matter, as well as its origin, have been investigated by several authors [|l9l -|25|l. In general, 
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one should expect that the short-range repulsion of the nuclear interaction originates spin 
alignment at sufficiently high densities; this has been verified in various model calculations 
extending from relativistic approaches p3J23] to HF description with finite range effective 



interactions |[25|| . The effect of a short range repulsion in neutron matter is most conveniently 
illustrated in the case of a Skyrme force without velocity dependent terms. In this case 
a transition takes place provided the density dependent repulsion grows rapidly enough. 
Indeed, a zero range interaction does not contribute to the potential energy in fully polarized 
neutron matter, because the Pauli principle forbids configurations of two neutrons with the 
same spin at the same point. As a result, when the repulsion in spin symmetric neutron 
matter exceeds (2 1 / 3 — 1) times the kinetic energy, the cost in kinetic energy requested 
to turn all spins up becomes balanced by the gain in potential energy, thus making the 
fully polarized state more favorable. The previous analysis also shows how to work out 
interactions without magnetic transitions. One needs a weak density dependent repulsion, 
more realistic to reproduce the monopole resonance, in order to reduce the gain in potential 
energy. To explore the results which would be obtained with such a force we have slightly 
modified the interaction SLy230b so as to increase the value of the critical density up to 
a factor 2, without impairing the equation of state. The results are indicated in Fig. 3 as 
SLy* for T = 5 MeV (solid line), T = 10 MeV (dashed line), and fcjC = 5 MeV in both cases. 
We find that the mean free path remains roughly constant up to densities around 0.9fm -3 . 

Further studies are obviously needed in the vicinity of the magnetic transition. Indeed, 
it occurs in a region where inputs from more realistic calculations are necessary for the 
determination of a reliable Skyrme type parametrization. Furthermore, other types of phase 
transitions, such as strangeness production or a quark-hadron transition, are likely to occur 
before one reaches this region. 
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FIGURES 





'"I V'l\ 

FIG. 2. Dynamical structure factors S^ s \oj,q) for S=0 and 1 channels (left and right parts, 
respectively), as functions of energy for interaction SIII at a temperature T = 5 MeV. Upper and 
lower parts respectively correspond to densities 0.1 fm" 3 and p c . Dashed and solid lines respectively 
correspond to transferred momenta q= 5 and 100 MeV/c. 



15 




0.1 0.2 0.3 0.4 0.5 0.6 

p(fm- 3 ) 



FIG. 3. The scattering mean free path for the neutrino as a function of neutron density for 
the parametrization SLy230b and its modificated version SLy*. Full and dashed lines respectively 
correspond to temperatures 5 and 10 Mev. In each group, from above to below the lines correspond 
to incoming momenta ki = 5, 10 and 3 T (in MeV/c). 
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FIG. 4. The relative opacity A~ 1 /Aj r 1 ee as a function of the density for interactions SIII and 
SLy230b. Full and dashed lines respectively correspond to T = h L c = 5 and 10 MeV. A horizontal 
line corresponding to relative opacity equal to unity has been drawn for easier comparison. 
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